Cherenkov radiation emitted by ultrafast laser pulses and the generation of coherent 

polaritons 



m 
o 
o 

(N 
C 

00 



o 

CZ3 



c3 



13 

o 
o 



> 

(N 

o 



o 
m 
o 

-I— > 



i 

C 

O 

o 



J. K. Wahlstrand and R. Merlin 
FOCUS Center and Department of Physics, The University of Michigan, Ann Arbor, MI 48109-1120 

(Dated: February 2, 2008) 

We report on the generation of coherent phonon polaritons in ZnTe, GaP and LiTaOs using ultra- 
fast optical pulses. These polaritons are coupled modes consisting of mostly far-infrared radiation 
and a small phonon component, which are excited through nonlinear optical processes involving 
the Raman and the second-order susceptibilities (difference frequency generation). We probe their 
associated hybrid vibrational-electric field, in the THz range, by electro-optic sampling methods. 
The measured field patterns agree very well with calculations for the field due to a distribution 
of dipoles that follows the shape and moves with the group velocity of the optical pulses. For a 
tightly focused pulse, the pattern is identical to that of classical Cherenkov radiation by a moving 
dipole. Results for other shapes and, in particular, for the planar and transient-grating geometries, 
are accounted for by a convolution of the Cherenkov field due to a point dipole with the function 
describing the slowly-varying intensity of the pulse. Hence, polariton fields resulting from pulses of 
arbitrary shape can be described quantitatively in terms of expressions for the Cherenkov radiation 
emitted by an extended source. Using the Cherenkov approach, we recover the phase-matching 
conditions that lead to the selection of specific polariton wavevectors in the planar and transient 
grating geometry as well as the Cherenkov angle itself. The formalism can be easily extended to 
media exhibiting dispersion in the THz range. Calculations and experimental data for point-like 
and planar sources reveal significant differences between the so-called superluminal and subluminal 
cases where the group velocity of the optical pulses is, respectively, above and below the highest 
phase velocity in the infrared. 

PACS numbers: 71.36,+c, 78.47.+P, 41.60.Bq 



I. INTRODUCTION 



Polaritons are hybrid modes of solids and, more gen- 
erally, dense systems associated with elementary exci- 
tations that carry electromagnetic polarization^^ Such 
excitations couple to light, resulting in a new mixed par- 
ticle called the polariton. The polariton description is, 
in some sense, an alternative to that involving the index 
of refraction; polaritons are dressed photons. Phonon 
polaritons, lattice vibrations coupled to infrared light, 
were first studied experimentally by spontaneous Raman 
scattering. 3 Soon after, mixing of coherent polaritons 
with visible laser light was reported for GaP4i The gen- 
eration, dispersion and decay of phonon polaritons have 
been extensively investigated using various coherent op- 
tical methods, particularly, coherent anti-Stokes Raman 
scattering (CARS) 5 followed by time-domain CARS 6 
and, more recently, the closely related technique of im- 
pulsive stimulated Raman scattering (ISRS). 

This work centers on ISRS, i.e., coherent scatter- 
ing using pulses whose bandwidth exceeds the polari- 
ton frequency. 7 In ISRS experiments, a pump pulse im- 
parts an impulsive force on the phonon component of 
the polariton starting a coherent oscillation which in 
turn, through the inverse process, perturbs the index of 
refractionA^ This perturbation is measured by a probe 
pulse that follows behind the pump pulse at a controllable 
time delay. The most common ISRS approach to gener- 
ate coherent polariton fields is to use an intensity grating 
obtained from crossing two pump beams. This geom- 



etry leads to polaritons of nearly well-defined wavevec- 
tor even in the pulsed case where the (transient) grat- 
ing travels with the pulsesii2*ii An alternative method, 
associated with the nonlinear susceptibility x^ 2 \ is to 
excite the electromagnetic component of the polariton 
through difference frequency generation (DFG), or opti- 
cal rectification. 12 Here, the pump pulse induces a non- 
linear polarization proportional to its intensity envelope. 
If the pulse is tightly focused and its group velocity is 
greater than the phase velocity in the infrared, its po- 
larization leads to emission of infrared light in much 
the same way that a relativistic dipole emits (coher- 
ent) Cherenkov radiation (CR). This interpretation, orig- 
inally proposed for second-harmonic generation^ was 
put forth for DFG by Auston and coworkersi^i 5 ^ 6 - whose 
technique became later the standard for generating short 
THz pulses. 17 . Recently, the CR interpretation was revis- 
ited in the analysis of pump-probe experiments on ZnSc 18 
and ZnTe^ In these and other polar materials, strong 
dispersion near the frequency of the transverse-optical 
(TO) phonon results in two qualitatively distinct regimes 
for CR depending on whether the velocity of the source is 
larger or smaller than the phase velocity of light at zero 
frequency. In the former, superluminal regime, the CR 
field is qualitatively the same as in the experiments of 
Auston and co-workersi& whereas the latter subluminal 
regime displays new features^ 

In the ultrafast literature, the relationship between the 
two, transient grating and tightly-focused single pump 
methods is rarely discussed. While the generation mech- 
anisms are closely related, ISRS is often described as 
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four-wave mixing and, as such, as a technique that dis- 
tinguishes itself from DFG. Here, we present a unified ap- 
proach encompassing the two methods. We show that the 
generation of THz radiation by single or multiple pulses 
of arbitrary shape can always be described as CR. We 
also report experimental results on three materials, ZnTe, 
GaP, and LiTa03, in three geometries (point, planar, and 
transient grating) using an electro-optic field sampling 
technique based on that of Auston and co- workers^ This 
method allows us to map out the polariton field as a func- 
tion of time delay and lateral position for polarization 
sources of arbitrary shape and, simultaneously, to image 
the shape of the source pulse itself. 



II. PHONON POLARITONS AND OPTICAL 
SUSCEPTIBILITIES 

As mentioned earlier, polaritons can be described ei- 
ther as excitations that carry polarization through the 
coupling to the electromagnetic field, or as light per- 
turbed by the refractive index of the medium, which ex- 
hibits strong dispersion near the resonant TO frequency. 
In the Cherenkov picture of phonon polariton generation 
we take the latter view, incorporating the effect of the 
infrared-active phonon into the dielectric function and 
the second-order nonlinear susceptibility. For a diatomic 
cubic lattice, and neglecting damping, the dielectric func- 
tion is of the Lorentz form 



e(fi) 



l-(fi/^T0) 2 ' 



(1) 



where J7to is the TO frequency of the three-fold degen- 
erate infrared-active mode, eoo is the dielectric constant 
including the effects of higher- lying resonances, and en is 
the static dielectric constant. The Lyddane-Sachs- Teller 
(LST) formula Q^o = 7 s " ^to gives the relationship be- 
tween r^LO, the longitudinal optical (LO) phonon fre- 
quency, and ^to- The refractive index n(fl) — y / e(TT), 
phase velocity c/n(f2), and polariton dispersion relation 
il = cq/n(£i) are plotted in Fig. ^ f° r a material with 
e(f2) given by Eq. In this work we deal mainly with 
frequencies below f^TO- Near zero frequency, where the 
dispersion curve is linear, the polariton is mostly light- 
like whereas, near Qto, it is phonon- like. When there 
is more than one infrared-active phonon, as in LiTa03, 
each mode contributes to the polarization and the dielec- 
tric function reflects a sum of these contributions (poles). 

Nonlinear interaction of visible light with phonon po- 
laritons involves both the lattice and the THz electro- 
magnetic field. The nonlinear tensor that couples po- 
laritons and visible light is thus composed of two parts: 
the Raman tensor ^ = x R , which couples visible light 
to the phonon mode with amplitude Q, and the second- 
order nonlinear susceptibility ^ = x^ > which couples 
visible to infrared electric fields. Both of these suscepti- 
bilities are almost entirely electronic in nature. Counting 




frequency " 

FIG. 1: Optical parameters for a generic material with a 
single triply degenerate infrared-active phonon; see Eq. Q 
(a) Refractive index n, (b) phase velocity c/n, (c) polariton 
dispersion. The darkened regions indicate the band between 
Qto and Q,lo in which light cannot propagate. 



both contributions, the total nonlinear polarization is 4 



P. 
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X ( ^ k E 3 {n)E k {u L ) + X f h QjE k {^) (2) 



where Qj is the amplitude of one of the three-fold degen- 
erate phonon modes, ldl is the frequency of the visible 
laser light and Q is the frequency of the phonon polari- 
ton. As is well-known, the Raman contribution can be 
included in an effective nonlinear susceptibility dijk, each 
tensor element of which is of the form 



o 2 , 

"TO 



(3) 



where the Faust-Henry coefficient C = 
X R (x^ 2 ^)~ 1 e*M~ 1 (^To)~ 2 gives the ratio of the Ra- 
man (often termed ionic) to the electronic contribution. 4 
Here, e* is the transverse effective charge and fj, is the 
reduced mass of the lattice mode. Because X K and x' 2 ' 
do not change appreciably with frequency in the THz 
range, C can be approximated as a constant for a given 
laser wavelength and material. If there is more than one 
phonon mode, there is a Raman contribution to d and 
to the Faust-Henry coefficient from each one. In our 
experiments we do not need to consider the dispersion 
of d, for we mainly deal with low frequency polaritons. 
Thus, in the following we assume that the nonlinear 
susceptibility is constant, and approximately given by 
X (2) (l + C7). 
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III. THEORY OF CHERENKOV EMISSION IN 
AN INFINITE CRYSTAL 

In non-centrosymmetric media, an optical pulse in- 
duces through a low-frequency polarization propor- 
tional to its intensity. This is DFG, also known as optical 
rectification^ The dipolar charge distribution moves at 
the optical pulse's group velocity v g and, like a charged 
particle, it can emit CR if v g is greater than the phase 
velocity of the THz radiation in the medium. The band- 
width of the THz radiation emitted by a pulse is limited 
to the pulse bandwidth. The spatial and temporal shape 
of the pulse envelope, together with the optical properties 
of the medium, determine the THz field. 

Kleinman and Auston adapted the theory of CR to 
an optical pulse traveling through a nonlinear isotropic 
medium^ Their approach, which we build on here, is es- 
sentially to separate the nonlinear mechanism by which 
the polarization in the medium is generated from the 
Cherenkov mechanism by which infrared light is emitted. 
This approach works well in the absence of higher-order 
effects such as self-focusing and self phase modulation, 
which affect the spatiotemporal shape of the pulse as it 
propagates, and depletion, which diminishes the induced 
polarization. Neglecting these effects, the emission pro- 
cess is entirely analogous to the emission of CR by a 
continuum of dipoles traveling at the group velocity of 
the pulse. In the following, we discuss electro-optic CR 
by pulses of a variety of spatiotemporal shapes. 

A. Point source 

For a beam focused to a small waist compared to the 
polariton wavelength, the polarization induced by the 
optical pulse can be approximated by a point dipole. 
This geometry, the 3D case discussed by Kleinman and 
Auston^ is the closest electro-optic analogue to conven- 
tional CR (emitted by a relativistic monopole) , which is 
coherent. The theory of CR was largely worked out more 
than 60 years ago by Tamm and FrankiSLSS For an in- 
finite medium and pulses propagating along the z axis, 
the fields are a function of z — v g t and the cylindrical co- 
ordinates p and (f>. This approximation, valid for a thick 
sample, greatly simplifies the theory. 

The radiation field due to a propagating dipole can 
be derived from the radiation field due to a monopole 
by taking a derivative in the direction of the dipole 
orientation^ In the experiments discussed here, the in- 
duced dipole is oriented perpendicular to the z axis. 
For this orientation, the fields are proportional to cos <fi, 
where (f> is measured from the direction of the dipole, 
dropping for simplicity the dependence, we haveSiiSi^ 

E z (p,z-v g t) <x fe^^i^— -1)|^M|4) 
J v z g n z (YL) dp 

E p {p,z-v g t) cx J ' -ie in{t -^ ] g^ dfl, (5) 
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FIG. 2: Calculated E p for a point dipole at p — and z — 
v g t — traveling through a medium with dispersion given by 
Eq. (a) Superluminal, v g > c/n(0), and (b) subluminal, 
v g < c/n(0), regimes. 



where E z and E p are the z and p components of the 
electric field, and 

a(p,Q) = H^[ S (Q)p] (6) 

where Hq 1 ^ is a Hankel function and s(f2) = Q(ev g /c 2 — 
If^/vg = Qtan9 c /v g . 

Radiation is emitted at frequencies for which the 
Cherenkov condition v g > c/n(f2) is satisfied. Inte- 
grating Eq. (J5J over these frequencies, two qualitatively 
different radiation patterns emerge, depending on v g m^ 
For v g > c/n(0), the laser pulse propagates faster than 
the entire lower polariton branch < ft < Uto- We 
call this the superluminal regime. For v g < c/n(0), the 
laser pulse only travels faster for frequencies in the band 
flc < Q < UtOi where the velocity matched frequency 
Q.C is defined by v g = c/n(flc)- We have referred to this 
as the subluminal regime^ 

Shown in Fig. [2] are plots of E p , calculated using 
Eq. |J5J, in a medium with refractive index given by 
Eq. (JIJ. For the superluminal regime (Fig. there 
is a sharp shock front with small amplitude ripples be- 
hind it. The weak ripples are due to dispersion, but the 
shock wave is defined by fi = 0. On the other hand, the 
subluminal regime (Fig. is dominated by dispersion 
effects. The angle of the phase front depends on p, and 
there is a large field near the axis of motion p = 0. The 
radiation pattern is still confined within a cone, but the 
cone angle becomes larger as v g decreases, in contrast to 
the superluminal regime. The special features of the sub- 
luminal regime are apparent in the calculations reported 
by Afanasiev and co-workers iSi 

In the CR literature, it is common to use the asymp- 
totic form of the fields for large p to find an approximate 
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frequency 



FIG. 3: Plot of angle versus polariton frequency; see Eq. 1101 
Dashed line indicates the Cherenkov angle 6c and fic indi- 
cates the group velocity matched frequency, (a) Subluminal 
regime, (b) Superluminal regime. 



analytical solution. In this limit, E„ becomes 



frequencies associated with this angle are solutions to 
s"(f2) = 0. In the superluminal regime, e(n)v 2 /c 2 — 1 > 
for all frequencies within the range of integration, and 
f2 = is the only solution. Thus, from Eq. llUfl . 
cos 8c — v g n(0)/c, which is the standard Cherenkov for- 
mula. In the subluminal regime, there are two solutions. 
For a material with dispersion given by Eq. (|T|). 9r — 
tan- 1 (4 7 -V('7-C)/C 3 ), where 7 = (1 - e^/c 2 )- 1 / 2 , 
C - 2 - (4 - Zr,) 1 / 2 , and 77 = 1 - (tt c /tt TO ) 2 3 

So far, we have considered isotropic materials, and 
therefore, the considerations above apply to ZnTe and 
GaP but, in the strictest sense, not to the uniaxial 
LiTaOa. The theory of CR in an anisotropic material 
is considerably more complex than that discussed here, 
because extraordinary and ordinary waves may be emit- 
ted. Other complications result for the case of an ex- 
tremely tightly focused pulse, for which there exists a 
small longitudinal component of the laser electric field, 
leading to the breakdown of the point dipole approxi- 
mation. In addition, the confocal parameter comes into 
play resulting in the breakdown of the infinite crystal ap- 
proximation. We note in passing that exact integrals for 
an electric dipole oriented arbitrarily with respect to the 
axis of motion, as well as those of a magnetic dipole, are 
given in the work of Zrelov^ 



E„ v f nT ° j^j e-iPNM-W^lffl, ( 7 ) 



O-c 



n 2 (n)^¥p 



For large p, the integral is small except for values of ft 
for which the integrand does not oscillate. Using the 
stationary phase method, Eq. 0can be approximated as 



i 



,-i[n i {t-z/v g )-s(n i )p] 



n 2 {n{)p y/\ s »(ni)\ 



where fii are the frequencies within the range of integra- 
tion for which the phase in the exponential is stationary. 
These are the solutions to 



ds 

'dn 



V a t — Z 



— cot I 



(9) 



which defines a cone of angle 9 with respect to the z 
axis. The cone angle depends on frequency because s(Q) 
contains n(f2). Explicitly, Eq. © becomes 



ds ^ 9 v a de v s 
dQ sc 2 ail il 



(10) 



The dependence of 9 on n is plotted for the two regimes 
in Fig. [3J In the superluminal case, 9 is a monotonic 
function of frequency, whereas in the subluminal regime 
two frequencies map onto the same angle, resulting in the 
complex beating behavior shown in Fig. [2p. 

The angle 9c , defined as the smallest angle at which ra- 
diation is emitted, is shown as a dashed line in Fig-EJ The 



B. Planar source 

When the beam waist w is much larger than the wave- 
length of the polariton (w » 2ncn/Q), we have another 
important limiting case. Here, the polarization induced 
by the optical pulse is shaped like a pancake. This planar 
geometry is often used for the generation and detection 
of THz pulses by nonresonant optical rectificationi2Si2i 
When applied to this problem, the CR picture is useful 
because it quantifies phase matching issues encountered 
in thick crystals due to dispersion in the infrared. 

For an infinitely thick crystal, and approximating the 
optical pulse as infinitely short (this is valid as long 
as the pulse is short compared to the period of the 
radiation being generated) the polarization induced is 
P ~ Q5{z — v g t)ed, where is a unit vector pointing 
in the direction of the dipole orientation and £ is the 
areal polarization. It is possible to obtain the CR field 
from the point dipole solution by convolution, but here 
we use the corresponding Hertz potential 



n(fi) = - n — e d . 



v g e(fl)s 2 (n) 
The electric field is theni& 
'' f2TO e(Q)n 2 



E 



-n(fi)dn 

4n( f QTO Sl 2 e ~ in ( t_2 / t '9) 



re,/ 



s 2 (fi) 



(ii) 

(12) 
(13) 
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FIG. 4: Calculated E p for an extended source of width w, ob- 
tained by convoluting results for a point dipole with a Gaus- 
sian in the p direction (subluminal regime), (a) w — 20 (jm, 
(b) w = 40 jLtm 
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FIG. 5: Calculated E p due to a transient grating, showing 
the checkerboard polariton field pattern generated behind the 
grating. The grating wavevector in (a) is smaller than in (b). 



It can be shown that the field vanishes unless s 2 (Sl) = 0. 
In the superluminal regime, s 2 (fl) > for all f2, and 
E = 0. In the subluminal regime, s 2 (ttc) = 0. Note that 
the phase matching condition is only met in the presence 
of dispersion. In the superluminal regime, no radiation is 
emitted from the bulk because the pulse is traveling faster 
than the phase velocities of all infrared frequencies. 

Practically, a beam cannot be focused to a point, nor 
is it really planar. If cylindrical lenses are used to fo- 
cus the beam, it can be extended in one direction and 
tightly focused in another. The extreme case of the 2D 
geometry was treated by Kleinman and Austen^ The 
point geometry solution is the Green's function of the 
problem, in that the radiation pattern of an arbitrary 
source can be obtained by convolution of the point solu- 
tion with the source shape. Two examples are shown in 
Fig. ^ This is a powerful machinery for calculating the 
radiation field generated by ultrafast optical pulses which 
have been shaped spatially and temporally, for example, 
in the manner recently demonstrated by Koehl, Adachi, 
and Nelson^ To the best of our knowledge, our approach 
has not been previously applied to this problem. 



C. Periodic source (transient grating) 

The last source shape we consider is a periodic dis- 
tribution, or transient grating, implemented by interfer- 
ing two pump beams inside the sample. This geometry 
is useful because it excites a unique polariton wavevec- 
tor determined by the grating period. If we cross two 
pump beams of equal intensity at an angle a, Gaussian 
in space with waists w x in the x direction and w y in the 
y direction and Gaussian in time with pulsewidth r, the 



intensity grating has the form 

I(x, y, z - v g t) oc e-^K^K cosf-^e"^-^ V- 2 

sin a 

(14) 

where kh is the wavevector of the laser light inside the 
medium. Here, we have assumed that a is small, which 
makes v g approximately equal to the group velocity of 
each separate pulse. This is a good approximation for the 
experiments described here, since k^ is much larger than 
the wavevector of the gratings. The Cherenkov angle con- 
verts the spatial periodicity of the transient grating into 
a temporal periodicity in the polariton field, resulting in 
a checkerboard pattern and a well-defined polariton fre- 
quency Q, as shown in Fig. [S] 

In the impulsive limit, a defocused probe pulse travel- 
ing at the same speed as the pump pulse and, thus, at the 
speed of the polariton field, samples a slice of the checker- 
board field in the xy plane. The electric field in this slice 
is periodic with the same wavevector k x = fc^/sina as 
the intensity grating. Through the electro-optic effect, 
the probe pulse diffracts from the periodically varying in- 
dex of refraction, resulting in an oscillatory signal at 2f2 
as a function of time delay, because the diffracted inten- 
sity is proportional to \E\ 2 . This technique is known as 
homodyne detection in the ISRS literature; 10 Heterodyne 
detection(iii22^S in which one mixes the diffracted signal 
with undiffracted probe light, gives a much larger signal 
at fl, because the signal is proportional to E. Because it 
is sensitive to the phase of the polariton field, which de- 
fines the phase of the diffracted probe pulse, heterodyne 
detection works just as well outside the checkerboard, in 
what we might call the Cherenkov wing of the field. As 
we discuss in detail below, heterodyne detection is often 
discussed in terms of four-wave mixing^ 

It is important to emphasize that the wavevector of 
the grating is not the same as the polariton wavevector, 
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since the polariton propagates in the direction normal to 
the wavefront. The magnitude of the polariton wavevec- 
tor is given by q = k x sin where 6 is the Cherenkov 
angle. In many early ISRS papers, it is assumed that the 
wavevector of the grating is the same as the polariton 
wavevector, tantamount to assuming that 6 is very large. 
This is generally a poor approximation, and leads to in- 
correct polariton dispersion measurements, as recently 
pointed out by Heblingi^i 

The simple interpretation above works best for large 
w, in which case the polariton wavevector is well-defined. 
For small w x , edge effects make the polariton held more 
complex but, by convolving the point dipole solution with 
the grating (as in Fig. EJ, we can calculate the polariton 
held due to a grating of arbitrary shape. 



IV. PHASE MATCHING 

One source of confusion when comparing the tran- 
sient grating, the tightly focused, and the planar geome- 
tries is that, at first glance, the phase matching con- 
ditions seem entirely different. Here, we clarify this 
point. An ultrafast pulse contains a continuum of fre- 
quencies and wavevectors, and the generated polariton 
wavevectors are determined by conservation of energy 
and momentum. The Cherenkov interpretation is, of 
course, consistent with phase matching. In all geome- 
tries, the nonlinear polarization contains terms of the 
form e «(k2-k 1 ).r- l ( W2 -^ 1 )^ wnere k 2 Wlj an( j U2 de- 
note wavevector and frequency components of the optical 
pulse. For the planar geometry, the wavevector compo- 
nents all point along the z axis, and conservation of mo- 
mentum requires q = fc 2 — k± = [w 2 n(a> 2 ) — win(wi)]/c, 
where q is the polariton wavevector. Since f2 = lu± — u>2, 

w 2 n(w 2 ) - (^2 + £l)n(w 2 + Q) , lr . 

q = (15) 

c 

sa — [n(u 2 ) + w 2 n (w 2 )] = — -. — r, (16) 

C V g {UJ2 ) 

making the connection to the group velocity matching 
argument derived in Section IIIB. 

A tightly focused optical pulse has wavevector com- 
ponents pointing in many directions. The dependence 
of the polarization on z — v g t forces q z = Cl/v g . The 
polariton wavevector and frequency must also satisfy 
Q = |q|c/n(n). Let q p be the component of the polari- 
ton wavevector perpendicular to the z axis. Momentum 
conservation gives 



'if. 



n/v g 



1, 



(17) 



and the polariton is emitted at an angle tan -1 q p /q z to 
the z axis, which is identical to the Cherenkov angle. 

In the transient grating geometry, the angle between 
ki and k 2 is well-defined, and the wavevector q = ki — k 2 



is emitted in the direction given by Eq. i|17|) because the 
polarization is a function of z — v g t. The transient grating 
configuration allows one to excite wavevectors that are 
higher than those accessible in the planar geometry. 

A recently proposed technique developed by Hebling 
and co-workers relies on a tilted pulse fronts This can 
also be encompassed in a Cherenkov radiation framework 
as a planar source which propagates at an angle. Tuning 
the tilt angle changes the wavevector of the generated 
polariton, and this shows promise as an easily tunable 
THz source. 

As is well-known, it is possible to interpret ISRS as 
four wave mixing^ Here one hides both the generation 
and detection of the polariton in xiffi much as we in- 
cluded the Raman contribution to the generation of the 
polariton in ■ The four-wave mixing polarization is 

(31 

P(ui 8 ) = ^(u p ,Q-u)p,u p )Ei(u)i)E 2 (wi -£l)E p (uj p ) 
= d(n)[d(Q)E 1 (u 1 )E^(u}i - n)]E p (u p ) (18) 

where E\ and £" 2 are the pump fields, E p is the probe 
field, and d is the effective nonlinear susceptibility defined 
in Section II. All propagation and damping effects can be 
included in \ c g , an d this formalism is often used£ 

Finally, we discuss briefly the advantages and disad- 
vantages of the various scattering geometries. The exper- 
iments with tightly focused beams produce the highest 
peak polariton amplitude (for a given pulse peak power), 
but also a wide range of wavevectors and frequencies. 
ISRS by a transient grating has the advantage of defin- 
ing the polariton wavevector and producing a large am- 
plitude polariton at that particular wavevector. CARS 
has the advantage of yielding the best frequency and 
wavevector resolution. 



V. MATERIALS 

Generation of electro-optic CR requires a non-zero 
X^ 2 \ present only in crystals without inversion symme- 
try. In our experiments, we used three readily available 
and well-characterized materials: the zincblende ZnTe 
and GaP, and LiTa03. The parameters relevant to our 
work are listed in Table [I] 



A. GaP and ZnTe 

Materials with the zincblende structure have a triply 
degenerate infrared-active optical mode. For light nor- 
mally incident on a (110) face, the selection rules forbid 
excitation of the nondispersive LO mode. For the TO 
mode, 
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(19) 
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Parameter GaP ZnTe LiTaQ 3 



«to/27t(THz) 11. 6 5.32 c 
n L o/27r(THz) 12. 1 6 6.18 c 
■JeE 3.3l f 3.16 c 



n a at 800 nm 3.56 / 3.24 9 
n c /27r' 1 (THz) 7.7 2.2 



12 e 
6.3 d 
2.2 



"Parameters are for the lowest A\ mode. 
b Reference 23 
c Reference M 
d Reference 3JJ 

c From LS T re lation using parameters of Ref. l35l 
^Reference [ 3(1 
9 Reference 1371 
^Calculated from £oo. 

TABLE I: Table of parameters used in the calculations. 



where R\ and i?2 are the Raman tensors for the modes 
polarized along the [110] and [001] directions, respec- 
tively. The single independent tensor component b corre- 
sponds to the nonlinear coefficient d^\. The dipolar dis- 
tribution induced by the pump pulse is always oriented 
perpendicular to the [110] direction. We measure the po- 
lariton by the inverse process, and we are only sensitive 
to the electric field component E p . 

In their pioneering work using nonlinear mixing of visi- 
ble and infrared cw laser light, Faust and Henry measured 
C = —0.47 in GaP. 4 This leads to total destructive inter- 
ference between the ionic and electronic contributions to 
d at 6.5 THz. Because this frequency is close to flc for 
Thsapphire pulses, this results in a small signal. Recently 
Leitenstorfer and co-workers found C = —0.07 in ZnTe 
by analyzing data from time-domain spectroscopy of THz 
pulses by electro-optic sampling^ This small value leads 
to d = at 5.2 THz, very close to J7to- For our experi- 
ment in ZnTe, we are only able to measure up to roughly 
4 THz, not high enough to probe effects of the dispersion 
of d. 



B. LiTaO; 



3 



This material crystallizes in the perovskite structure. 
LiTaC>3 is ferroelectric at room temperature and exhibits 
4 A\ symmetry and 9 E symmetry modes, all of which are 
Raman and infrared active^ There has been consider- 
able disagreement over the assignment of several Raman 
peaks, but the lowest-frequency A\ mode considered here 
is well-characterized. The Raman tensor for A\ phonons 

i8 aa 




(20) 



The independent tensor elements a and b are associ- 
ated with nonlinear coefficients d 3 \ and c?33, respectively. 



Here, we use \^ — d 33 and x R = b, which are the 
largest coefficients, to both generate and detect the po- 
lariton. The dipolar distribution points along the optic 
axis, and the probe pulse is sensitive to the component of 
the Cherenkov electric field along that axis. Barker, Ball- 
man, and Ditzenberger measured the infrared reflectivity 
and found that, for the extraordinary index, most of the 
oscillator strength in the infrared is in the lowest-lying 
A\ TO mode^ This mode is also the strongest Raman 
scatterer. We performed our experiment on the lowest 
polariton branch. 

The refractive index in the visible is weakly 
anisotropic. Since the laser pulses in our experiment are 
polarized along the optic axis, there is no effect of this 
anisotropy, in contrast with experiments on E modest 
The static dielectric constant is also weakly anisotropic. 
However, since the lowest-lying A\ and E TO frequen- 
cies are different (6 and 4.2 THz, respectively), the di- 
electric function becomes more anisotropic at large fre- 
quencies, where the cubic approximation is expected to 
break down. For the low frequencies discussed here, the 
approximation is well obeyed. For the lowest-lying A\ 
mode in LiTaOs, C is positive^ 



VI. EXPERIMENTS 
A. Techniques 

Our experimental technique is similar to that of Aus- 
ton and Nuss^ A diagram of the experimental setup is 
shown in Fig. EJ Pump and probe beams are focused 
and crossed inside the sample. The pump pulse is the 
Cherenkov source and the tightly focused probe serves 
as the detector of the Cherenkov radiation field. The 
probe pulse "surfs" on top of the Cherenkov wakefield 
generated by the pump pulse and samples the polari- 
ton field through both the linear electro-optic and the 
Raman effects, which both cause changes in the refrac- 
tive index. We change the time delay (t) between pump 
and probe with a motorized delay stage in the probe 
path and the relative focal position (p) between pump 
and probe by moving a motorized translation stage upon 
which the pump's focusing lens is mounted. By sampling 
the diffracted probe beam as a function of t and p, we 
map the polariton field. 

We used two different detection schemes. The first is 
a frequency shift measurement, which takes advantage of 
the fact that the spectrum of the probe is shifted pro- 
portional to dn/dt& The polariton induces a change in 
the index of refraction due to its electromagnetic and 
vibrational components, so that dn/dt contains a term 
proportional to dE/dt and another one to dQ/dt. We 
measure the frequency shift by cutting half of the probe 
spectrum with a bandpass filter before detection with a 
photodiode. This technique is especially useful when the 
pump and probe polarizations are chosen perpendicular 
to each other to suppress the scattered pump light. The 
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FIG. 6: Experimental setup. Pulses split into pump and probe beams, the time delay t between the pulses is varied with an 
automated translation stage. The relative position of the two beams p is adjusted by moving a single lens. The inset shows a 
close-up of the geometry inside the crystal. The pump and probe pulses travel along parallel paths inside the crystal at the 
same group velocity v g . BS and PBS denote, respectively, a 10:1 beam splitter and a polarizing beam splitter. Not shown are 
A/2 waveplates, used to rotate the polarization of the pump and probe beams, and a A/4 waveplate in the probe path which 
resolves the induced birefringence due to the electro-optic effect. 



second technique, widely used for detecting THz radia- 
tion, is polarization detectioniSSi^ When properly ori- 
ented, the polariton's electric field induces birefringence, 
which can be probed by measuring the polarization state 
of the probe beam after the sample. This yields a signal 
directly proportional to the THz Chcrenkov electric field. 

The spot size inside the crystal was determined ex- 
perimentally using the same technique. In ZnTe and 
GaP, two photon absorption depletes the probe when 
pump and probe overlap temporally and spatially, creat- 
ing a large Gaussian peak in the signal. In LiTa03, the 
intensity-dependent refractive index due to x*- 3 ^ modu- 
lates the probe pulse's phase when it overlaps with the 
pump pulse, also producing a large signal (known as the 
coherent artifact). 

Auston and Nuss used thin samples and a collinear ge- 
ometry to achieve spatial resolution on the order of 1.5 
/j,mm& A problem with this setup is that it is difficult to 
measure the field behind the source, because of strong 
pump scattering due to the fact that the two beams are 
collinear. In order to measure the source itself plus the 
field behind it, we used the nearly-collinear setup de- 
scribed above, which enables the spatial separation of the 
pump and probe beams even when they overlap spatially 
inside the sample. This introduces a source of distortion 
in the data due to spatial walk-off between the pump 
and probe pulses as they travel through the sample. To 
minimize this walk-off, we focused the pump and probe 
beams with two small lenses placed as close together as 
possible. In order to ensure that the walk-off did not 



affect the spatial resolution, we made the plane shared 
by the pump and probe beams perpendicular to the axis 
along which the pump lens was moved. We accomplished 
this by directing the probe beam below the pump lens, 
and moving the relative focal position horizontally. The 
resolution in z — vt is determined by the pulsewidth of 
the probe pulse, and in p by the focal waist of the probe 
beam. 



Ultrafast pulses were generated by a Tksapphire oscil- 
lator producing 800 nm, 60 fs pulses at the repetition rate 
of 82 MHz. Acousto-optic amplitude modulators were 
used to modulate the pump and probe beams at 3 MHz 
and 2.95 MHz, respectively, and the probe was detected 
with an unbiased photodiode at the difference frequency, 
50 kHz, with a lock-in amplifier. The maximum pump 
and probe average power at the sample were ~140 and 
30 mW, respectively, limited mostly by the efficiency of 
the modulators. We emphasize the importance of the 
double-modulation technique because it helps reject the 
scattered pump light. Recently, experiments have been 
performed using a Zernike method with a CCD camera, 
to image the phase change of a weakly-focused probe 
beam caused by the polariton field as a function of po- 
sition and time delayiS&iii This is another powerful tech- 
nique for imaging the Cherenkov field, although it suffers 
from unwanted large scattering near the pump pulse. 
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FIG. 7: (a) Experimental results for LiTaOa. The pump waist 
is 18 pm. (b) Convolution of the point dipole field with a 18 
/im Gaussian. 

B. Point source: Experiments 

If the pump pulse is focused to a waist smaller than 
the wavelength of the polariton, we can approximate it 
as a point dipole. Because of the large thickness of our 
samples, we are limited to a focal waist of about 20 /im 
which corresponds to a maximum polariton frequency of 
2-4 THz, depending on the material. For lower frequen- 
cies, the pump pulse acts like a point dipole. 

First we discuss the data on LiTaOs, a superluminal 
material; see Fig. \7\ Our sample was an x-cut wafer 2 
inches in diameter and 1 mm thick, and we used fre- 
quency shift detection. For pulses of wavelength 800 nm, 
Be = 69° . The experiment matches the calculations quite 
well after taking into account the 20 /im pump waist by 
convoluting the point dipole solution with a Gaussian of 
that width. However, the absence of ripples behind the 
sharp Cherenkov cone indicates that damping, not in- 
cluded in the simulation, is strong for higher frequencies. 
In addition to polaritons of the lowest branch, we were 
also able to generate and detect polaritons and coherent 
phonons from higher branches; see later. 

Our measurements on ZnTe, a subluminal material, 
and the comparison with calculations are shown in Fig. [8] 
The (110) oriented ZnTe sample was 10x10x1 mm 3 , 
and polarization-sensitive detection was used. The data 
shows many of the features predicted for subluminal CR. 
The artifact that appears about 1.6 ps after the pump 
pulse is due to a second pulse about 1000 times smaller 
than the pump pulse, due to a reflection from an optical 
element upstream. We have subracted the contribution 
from an exponentially decaying signal caused by free car- 
riers generated by two-photon absorption. This signal is 
partially balanced out by using polarization-sensitive de- 
tection, and the remnant has been fitted and subtracted 
out in Fig. [SJi to show the absence of CR near p = 0. 

As discussed earlier, subluminal CR is expected to ex- 
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FIG. 8: (a) Results of a pump-probe differential transmission 
experiment in ZnTe. An overall constant background was 
subtracted for each time scan, and a decaying exponential 
was subtracted out for time scans near p — 0. The feature 
1.6 ps behind the pump pulse is an artifact, (b) Calculation 
for a point dipole convoluted with a 22 pm Gaussian in the 
transverse direction, (c) Experimental data for a 300 pm thick 
ZnTe sample. 

hibit a beating pattern. In Fig. there is a striking 
node that separates the shock wave into two and possi- 
bly three distinct parts. This can be explained in the 
calculations by both subluminal beats, or by a convolu- 
tion artifact due to the finite spatial size of the pump 
and probe pulses. However, note that the node is not as 
prominent in the calculation. This discrepancy, and the 
lack of experimental signal near p = 0, is probably due to 
the finite size of the crystal, as discussed in detail later. 

C. Planar source: Experiments 

For the planar geometry, scans were performed only 
along the p axis in the z — v g t direction. Results for 
GaP and ZnTe with a pump waist of roughly 100 /im are 
shown in Fig. [5] The GaP sample, (110) oriented, was 
a 4x1x1 mm 3 crystal. The time scans show oscillations 
at the group velocity matched frequencies predicted by 
theory, with an apparent decay probably due to a com- 
bination of the finite length of the sample (see below), 
and the finite extent of the pump pulse, — The planar ge- 
ometry is commonly used for generating THz pulses by 
nonresonant optical rectification^ The distortion of the 
electro-optic sampling signal due to infrared dispersion, 
especially in a thick crystal, has been the subject of much 
discussion recentlyi^^ 

We also performed an experiment on the superluminal 
material LiTaOs. As expected, there was no signal from 
the lowest polariton branch. However, our pulses were 
short enough to see higher branches. The phase matching 
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FIG. 9: Results of pump-probe experiments with weakly fo- 
cused ultrafast pulses in (a) ZnTe and (b) GaP. The region 
near time zero is dominated by two photon absorption, (a) 
1 mm thick ZnTe crystal. The main frequency is roughly 2 
THz. (b) 1 mm thick GaP crystal. The main frequency is 
roughly 7.5 THz. The feature at 2.2 ps is an artifact due to 
a reflected pulse. 



condition is given by Q/v g , and this expression was found 
to match the higher polariton branches quite well. 



D. Finite crystal effects 

Point geometry experiments were also performed on a 
300 /im thick ZnTe sample. As is apparent in Fig. QJ, 
the node we observed in the 1-mm-thick crystal has now 
disappeared. This is probably a result of the breakdown 
of the infinite crystal approximation. In the superluminal 
regime, the infinite crystal calculation works well because 
the radiation is emitted at large angles at all frequencies 
down to zero. Therefore, effects of the finiteness of the 
crystal can be expected only at the very edge of the shock 
wave, at p = L tan (9(7 = Luq/c, where L is the length of 
the crystal. For example, the edge of the shock wave is 
at p 2.6 mm in our 1 mm thick LiTaC>3 sample, well 
beyond what we can practically measure due to the fall- 
off proportional to p -1 / 2 in the electric field from a point 
source. 

The situation is considerably different in the sublumi- 
nal regime because polaritons near £lc are emitted in the 
near forward direction. In the 1 mm thick ZnTe crystal, 
the Cherenkov angle for 2.5 THz polaritons is about 9°, 
so that this component reaches p — 160 /im before the 
pulse reaches the end of the crystal. In the 300 /im thick 
crystal, it only reaches p = 50 /im. Thus, the subluminal 
beats, due to interference between frequency components 
sent in the forward direction, require some propagation 



time to develop. We believe this is the reason why they 
do not appear in the thin crystal. This could also ex- 
plain other discrepancies between theory and experiment 
in Fig. El The large oscillatory signal at tie near p = in 
the calculation is clearly absent in the experiment, prob- 
ably because radiation at Qc is emitted in the forward 
direction The difference in the shape of the cone between 
theory and experiment is likely due to the same effect. A 
finite crystal calculation is probably needed to fully un- 
derstand the details of the experiment. Note that CR in a 
finite crystal is similar to the so-called "Tamm problem" 
of CR for a particle traveling along a finite path—ilii 

The finiteness of the crystal is expected to play an 
even more important role in the planar geometry, in 
which Cherenkov components are emitted exclusively in 
the forward direction. The CR pattern is the result of 
interference among these components, so that this pat- 
tern should change dramatically with crystal thickness. 
A very thin crystal would show no propagation effects, 
simply reproducing the time derivative of the pump pulse 
envelope as a short THz pulse. 



E. Periodic source: Experiments 

We adapted the setup described previously to perform 
experiments with two pump pulses, creating a periodic 
dipolar distribution in LiTa03, and used frequency shift 
detection to map the polariton field. The pump pulse 
was split into two equally intense pulses which were fo- 
cused by a single lens and crossed inside the sample at an 
angle a. Using a beam waist such that w << v g r/ sina, 
we avoided the so-called "pancake effect,"— a shrinking 
of the overlap between weakly-focused pulses due to their 
short duration. In our experiments, we did not have to 
measure a, because the grating period was directly ob- 
tained by measuring the coherent artifact as a function of 
spatial position. As for the point geometry, we mapped 
the polariton field due to this pump configuration as a 
function of z — v g t and p. The results in Fig. show 
excellent agreement with the theoretical predictions. 

F. LiTaOa: Phonons and Polaritons 

LiTaOs has a particularly rich history in the ISRS 
literature, 10 . 11 ! 28 ! 30 ! 43 . 47 ! 48 ! 49 ! 50 ! 51 ! 52 In this section we 
compare our data to those of other reports. There are 
two main points of contention in the literature: the ori- 
gin of frequency-dependent damping effects and avoided 
crossings at low frequency. Our results for the polariton 
dispersion are summarized in Fig. 1111 

As previously discussed for the transient grating geom- 
etry, see below Eq. ((T7|l . the Cherenkov angle appears in 
the expression for the polariton wavevector. Using the 
correct expression, our data fall almost exactly on the 
dispersion relation of the lowest branch A\ polaritons in- 
ferred from the IR data of Barker et al. 35 Recent ISRS 
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FIG. 10: Results of transient grating experiments in LiTa03. Left: Derivative of the polariton field as a function of time delay 
and position of the pump lens. Right: Cross sections at zero time delay, showing the intensity grating's period. The phase 
matched frequency is (a) 1 THz (b) 2 THz (c) 3 THz. Damping is higher at higher polariton frequencies. 



papers have begun to take this into account as welli22i2i 
The grating technique is also sensitive to non-dispersive 
phonons, and the signal at 7.8 THz roughly matches an 
A\ phonon peak seen in Raman spectral Scans in the 
planar geometry, with wavevector given by the group 
velocity of the pulse, are also shown in Fig. ^2 Here, 
we are sensitive to even more phonon branches, includ- 
ing a higher-branch polariton above 12 THz. Note that 
higher branches were also observed in a recent hetero- 
dyne experiment^ 

Our measurements of the polariton damping rate 7 
by the grating method^ agree very well with those re- 
ported by others: 7 = 2.4 cm -1 at 1 THz, 7 = 4.4 
cm -1 at 2 THz, and 7 = 7.5 cm -1 at 3 THz. Large 
damping, peaked at 2.7 THz, was reported by Auston 
and Nuss in their original experiment. They suggested 
that this was due to multi-phonon absorption. In their 
homodyne transient grating experiment, Wiederrecht et 
aU£ measured increased damping peaked at 2.76 THz and 
4.38 THz, and attributed this to strain-induced coupling 
to heavily damped iS-symmetry modes at the same fre- 
quencies. Very recently, Crimmins, Stoyanov, and Nelson 



measured the damping rate using heterodyne detection 
and suggested that the damping is due to quadratic cou- 
pling between the phonon-polariton mode at 2.1 THz and 
the lowest-lying .E-mode at 4.2 THz^S 

The inelastic neutron scattering data of Cheng et al 5 ^ 
supports the idea that the damping is due to anharmonic- 
ity, more precisely, to a process involving the simultane- 
ous creation and destruction of two different phonons. 
Close to the observed peaks at 2.76 and 4.38 THz, the 
phonon dispersion shows TO and LA (TA) branches for 
which the frequency difference at the zone edge, where 
the density of states is large, is 2.6 THz (4 THz). In 
this context, we note that decay by two-phonon dif- 
ference modes has been seen in many materials, in- 
cluding ZnTeS Furthermore, many spontaneous Raman 
studies22^&, including one on our particular crystal, 5 '' 
show a broad peak at 2.55 THz, consistent with a second- 
order feature. An alternative explanation for polari- 
ton broadening, supported by the fact that the TA and 
LA branches at the zone edge are at 3.2 and 4.6 THz, 
respectively, 54 is coupling to high-wavevector phonons 
due to disorder-induced breakdown of momentum conser- 
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FIG. 11: Summary of results in LiTaOs. Circles indicate fre- 
quencies measured using the grating technique, with wavevec- 
tor given by k x / sm6, where k x is the wavevector of the grat- 
ing and 9 is the Cherenkov angle. Crosses indicate measure- 
ments with a planar geometry, with wavevector given by Q/v g . 
The solid lines give the polariton dispersion according to the 
infrared measurements of Ref. [3^. 

vation. Given that LiTa03 is often non-stoichiometric, 
defect-related processes may contribute strongly to the 
damping. 

Another anomaly reported in LiTa03 is an extra 
avoided crossing in the polariton dispersion at 0.9 
rpjj Z( ii j 47 i 48 Similar avoided crossings were also reported 
for LiNb03j£& and increased absorption was observed at 
the corresponding frequencies in stimulated Raman gain 
measurements^ The avoided crossings were assigned by 
Bakker and co-workersii to quantum beats between lev- 
els of an intrinsic anharmonic potential. However, our 
work on LiTaC>3, as well as a recent heterodyne ISRS 
experiment,-- show no sign of an avoided crossing at 0.9 
THz. The observation of several low-lying modes in a 
CARS study of MgO-doped LiNbO^ suggests that the 
extra modes are due to defects. If this were the case, the 
reported discrepancies for LiTaOs may simply be due to 
differences in the concentration of defects for different 
samples. 

We take issue with the intrinsic anharmonicity pro- 
posed by Bakker et alii to explain these extra modes. 
Unlike in molecules, lattice anharmonicity in solids leads 
to decay, because anharmonic coupling involves a contin- 
uum of modes. Hence, unless the anharmonicity is strong 
enough to support 2-phonon bound states, it cannot in- 
troduce new frequencies. To the best of our knowledge, 



no such bound states have been observed in LiTaOs or 
LiNb0 3 . 

An unexplored area is the role of anisotropy in LiTaOs. 
As briefly discussed earlier, for frequencies approaching 
the lowest iJ-phonon frequency, the index of refraction 
becomes extremely anisotropic. Anisotropy is known to 
have an important effect on CR 23 but as far as we know, 
it has been completely ignored in the ISRS literature. For 
a point charge traveling through a uniaxial material in a 
direction perpendicular to the optical axis, both ordinary 
and extraordinary waves are emitted, with the extraordi- 
nary waves polarized at a finite angle to the optical axis. 
Assuming this is also true for a dipole oriented along the 
optic axis, it should have a large effect on the CR pattern. 

Finally, we note the observation of wavevector over- 
tones in transient grating experiments using amplified 
pulses, which Nelson et al^2S ascribed to lattice anhar- 
monicity. Since polariton nonlinearity is primarily due 
to phonon anharmonicity and the polaritons in question 
are mainly light-like, we consider this interpretation to be 
questionable. Instead, we propose that the higher-order 
diffraction peaks simply reflect a non-sinusoidal grating 
due to higher-order optical nonlinearities. This is an area 
which needs to be explored further. 



VII. CONCLUSIONS 

In summary, we have used ultrafast optical pulses to 
generate and image coherent polariton fields in ZnTe, 
GaP and LiTaOs. For single and multiple pulses of arbi- 
trary shape and, in particular, for the transient-grating 
geometry, we have shown that the polariton field can be 
calculated by convolving the slowly-varying pulse enve- 
lope with the Cherenkov field due to a dipole. The stan- 
dard phase matching arguments can also be explained by 
applying ideas from the theory of Cherenkov radiation. 
Results for point-like and planar sources reveal impor- 
tant differences between the superluminal and sublumi- 
nal regimes, especially for the Cherenkov emission by a 
plane of dipoles which can only occur in the subluminal 
case. 
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